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O ; Measurements of ep scattering with squared 4-momentuni transfer Q^ up to 35000 

D . GeV^ are compared with the expectation of the standard deep-inelastic model of 

'^ ' lepton-nucleon scattering (DIS). For Q'^ > 15000 GeV^, Nobs = 12 neutral current 

candidate events are observed where the expectation is N^js = 4.71 it 0.76 events. 
^ . In the same Q^ range, N^hs = 4 charged current candidates are observed where 

H I the expectation is Ndis = 1.77 it 0.87 events. The probability T'{N > Nobs) that 

the DIS model signal N fluctuates to A^ > Nobs in a random set of experiments 
is 6 X 10^^ for neutral current and 0.14 for charged current. The difference in the 
observed and expected number of Neutral Current events is mostly due to events at 
large masses M = y/xs in which the positron is backscattered at large y = Q'^ /M"^ . 
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1 Introduction 

The ep collider HERA offers the unique possibility to probe the proton at very small 
distances (~ 10~^^ cm) via t-channel exchange of highly virtual gauge bosons |l|, |^, 0, 
and to search in the s-channel for new particles which couple to lepton-parton pairs with 
masses up to the kinematic limit of i/s ~ 300 GeV 0]. 

In this paper, measurements of neutral current (NC) and charged current (CC) deep- 
inelastic scattering (DIS) are considered using all available e^p data collected by HI at 
HERA from 1994 to 1996. The positron beam energy E^ was 27.5 GeV and the proton 
beam energy Ep was 820 GeV. The total integrated luminosity amounts to 14.19 ± 0.32 
pb~^, an increase of a factor ~ 5 compared to the above cited earlier studies. 

2 The HI Detector 

A detailed description of the HI detector can be found in |Q . Here we describe only the 
components relevant for the present analysis in which the final state of the events involves 
either a positronQ with high transverse energy or a large amount of hadronic transverse 
energy flow. 

The positron energy and angle are measured in a liquid argon (LAr) sampling calorime- 
ter 1^ covering the polar angular^ range 4° < 6 < 153° and all azimuthal angles. The 
granularity is optimized to provide fine and approximately uniform segmentation in labo- 
ratory pseudorapidity and azimuthal angle 0. It consists of a lead/argon electromagnetic 
section followed by a stainless steel/argon hadronic section. Electromagnetic shower ener- 



gies are measured with a resolution of a{E)/E ~ 12%/ J E/ GeV© 1% and pion induced 



hadronic energies with a{E)/E ~ 50%/ J E/ GeV © 2% after software energy weighting. 
These resolutions were measured in test beams with electron energies up to 166 GeV ||^, |^ 
and pion energies up to 205 GeV p. The absolute energy scales are known to 3% and 
4% for electromagnetic and hadronic energies respectively. The angular resolution on 
the positron measured from the electromagnetic shower in the calorimeter varies from 
~ 2 mrad below 30° to < 5 mrad at larger angles. A lead/scintillating-fibre backward 
calorimeter |p extends the coverage^ at larger angles (155° < 6 < 178°). 

Located inside the calorimeters is the tracking system which is used here to deter- 
mine the interaction vertex. The main components of this system are central drift and 
proportional chambers (25° < 9 < 155°), a forward track detector (7° < 6* < 25°) and a 
backward drift chamber^. The tracking chambers and calorimeters are surrounded by a 
superconducting solenoid providing a uniform field of 1.15 T parallel to the z axis within 
the detector volume. The instrumented iron return yoke surrounding this solenoid is 
used to measure leakage of hadronic showers and to recognize muons. The luminosity is 
determined from the rate of the Bethe-Heitler ep -^ epy bremsstrahlung measured in a 



"'^The analysis does not distinguish explicitly between e~^ and e^. 

^The z axis is taken to be in the direction of the incident proton, and the origin of coordinates is the 
nominal ep interaction point. 

■^The detectors in the backward region were upgraded in 1995 by the replacement of the 
lead/scintillator tile calorimeter |1^ and a proportional chamber. 



luminosity monitor. This consists of a positron tagger (e-tagger) and a photon tagger 
(7-tagger), located —33 m and —103 m respectively from the interaction point. 

3 Event Selection and Kinematics 

3.1 Selection of ep Collisions 

The analysis considers the accumulated data corresponding to HI running conditions 
for which the central jet chambers (CJC), the LAr calorimeter and associated triggers, 
the backward calorimeter, the time-of-flight system and the luminosity system are fully 
operational. The trigger requirement is either an electromagnetic cluster or a large imbal- 
ance in transverse "momentum" -P/^iss measured with coarse trigger towers of the LAr 
calorimeter 0. Every tower k provides a fast energy sum Ek which contributes to the 
reconstruction of PrZiss = \J{T.E^^kY + {T^Ey^kY. 

Background not related to e^p collisions is first rejected by requiring for each event : 

1. a primary interaction vertex in the range | z — z \< 35 cm where z varies within 
±5 cm around 2; = depending on the HERA beam settings; 

2. that the event time to determined accurately (o"'""^'"(to) — 1-6 ns) from the central 
jet chambers coincides with the nominal time of the bunch crossings of the beams; 

3. that it should survive a set of halo and cosmic muon filters 0. 

Cut (1) mainly suppresses beam-wall and beam-residual gas interactions. Cuts (2) and 
(3) eliminate any remaining background produced by cosmic rays and by "halo" muons 
associated with the proton beam. 

3.2 Neutral Current Selection and Kinematics 

The following criteria are designed to select NC DIS events in the high Q^ domain. The 
cuts, which rely only on calorimetric information, require : 

1. an isolated positron with Ex^e > 25 GeV (-Et.e = EgSinOe), found within the polar 
angular range 10° < 9e < 145°; the positron energy cluster should contain more than 
?o of the LAr energy found within a pseudorapidity-azimuthal cone of opening 



(A?7e)2 + (A0e)^ = 0.25 where 77^ = -Intanf ; 
a total transverse momentum balance PT,miss/ J E^^e < 3vGeV where Px^r 



(S Ex^i) + iJ2 Ey^i) summed over all energy deposits i in the calorimeters {E^^i = 
Ei sin 9i cos 0j and Ey^i = Ei sin 9i sin 0j) ; 

3. a limited momentum loss in the direction of the incident positron, 43 <J2{E — Pz) < 
63 GeV, where E{E- Pz)=i: (Ei - Ez,i) with Ez,i = Ei cos Oi. 



The identification of positron induced sliowers relies on detailed knowledge of the expected 
lateral and longitudinal shower properties 0. The efficiency for the detection of positrons 
exceeds 90% everywhere within the acceptance cuts, the main losses being due to showers 
developing through the inactive material between calorimeter modules. The cut (2) makes 
possible a very efficient NC DIS selection up to the highest Q^ by taking into account the 
natural scale of the hadronic energy resolution. The cut (3) retains 90% of NC DIS events 
and exploits the fact that by energy-momentum conservation, the !](£' — P^) distribution 
for DIS events is peaked at 2ii^g. It rejects events where a very hard collinear 7 is emitted 
by the initial state positron. 

To calculate the appropriate DIS Lorentz invariants Q^, y and M = -/xs, two different 
sets of estimators are used. Firstly, use is made of only the measurement of the scattered 
positron: 

,. ^ /Ql" „2 ^T,e . E.-E.COS^. 



ye' '^ 2^0 

This method will henceforth be denoted the e-method. Secondly use is made of the 
reconstructed angles 9e and 9h of the positron and hadronic final state. Using ae = 



tan(^,/2) = {E,-E,^,)/ET,e: and a^ = tan(^j2) = E(^i -^.,*)/V(E^x,i)' + (E^y,^)^ 
where the summations are over all energy deposits of the hadronic final state, then [^ : 



M, 



2a 



E'e ^2 El 



Qia = -lt—< \ T' 2/20 



\aeahEp Epae{ae + ah) cxe + ah 



This method will henceforth be denoted the 2a-method. In the parton model, the variable 
X is interpreted as the momentum fraction carried by a parton in the proton. The centre 
of mass energy of the positron-parton collision is M. The variable y is related to the polar 
scattering angle 9^ of the positron relative to the incident proton in the centre of mass 
frame of the positron-parton collision according to y = (1 + cos9*)/2. 

In addition to the cuts (1) to (3) above, the analysis is restricted to the kinematic 
range Q^ > 2500 GeV^ and 0.1 < y < 0.9. The resolution in Me degrades with decreasing 
ye {5Mf./Mf. oc l/^/e) and so the low y domain is excluded. Excluding the high y values 
suppresses the photoproduction background (see section |^) and also avoids the region 
where QED radiative effects are largest for the e-method. For this kinematic range, the 
NC trigger efficiency exceeds 96% and is consistent with 100% to within experimental 
error. 

Following the NC selection described above, 471 (493) DIS event candidates are ac- 
cepted using the e-method (2Q;-method) for the calculation of the kinematic cuts. 

3.3 Charged Current Selection and Kinematics 

The selection of event candidates for the CC DIS sample requires: 
1. the total missing transverse momentum Pr.miss = PT,h > 50 GeV; 



2. the total transverse energy Et = Y. JEl^ + Ely calculated from energy deposits in 
the calorimeters to match P^.h such that {Et — PT,h)l Et < 0.5. 



These cuts eliminate the photoproduction and NC DIS background. 

The DIS Lorentz invariants Q^, y and M are calculated using the Jacquet-Blondel 
pj]] by summing over all measured final state hadronic energy deposits using: 



ansatz 



Mh 
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Vh ' 



Ql 
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EiE- P,) 



2E0 



This method will henceforth be denoted the /i-method. The /i-method can also be used 
for NC DIS. 

In addition to the cuts (1) and (2) above, the analysis is restricted to the kinematic 
range i/h < 0.9. The resolutions in both Mh and Ql degrade with increasing y (both 
SMh/Mh and 5Q\/Q\ behave as oc 1/(1 — yh) for yh ~ 1) and so the high yh domain is 
excluded. The requirement of cut (1) restricts implicitly the kinematic range to Q\ > 
2500 GeV^. Throughout this range, the CC trigger efficiency is > 90%. 

The CC selection described above retains 31 CC DIS event candidates. 



4 Standard Model of Deep-Inelastic Scattering 

The calculation of an expectation for NC and CC DIS ep scattering is performed using 
the parton model in the approximation of single 7/Z and W boson exchange. Detailed 
expressions for the NC and CC cross-sections which contain a description of the proton 
in terms of scale dependent structure functions can be found in 0]. For NC DIS, the con- 
tribution of the longitudinal structure function Fl is neglected. The structure functions 
are expressed in terms of parton densities which are taken here from the Mart in- Roberts- 
Stirling MRS (H) |T^ parametrizations. These were determined by a global fit of struc- 
ture functions measurements in fixed target experiments 
production and direct photon production 



14 



in 



p!7| , inclusive lepton 
and low Q^ structure functions 
measurements at HERA |23, |2^ . The parton densities are evolved to the high Q^ relevant 



18, 19, PO 2 



for this analysis using the next-to-leading order DGLAP equations p4[ and divergences 
are regulated in the DIS renormalization scheme [^]. Higher twist contributions are 
neglected. 

For the comparison with data, a Monte Carlo simulation following such an approach is 
performed using the event generator DJANGO [EH]. This generator includes the QED first 



order radiative corrections, and emission of real bremsstrahlung photons |2^ . Moreover 
it offers an interface to two possible models for the emission of QCD radiation. The 
LEPTO [ES] model includes the QCD matrix elements to first order in as, complemented 
by leading-log parton showers to model higher order radiation. The ARIADNE [^ 

to simulate QCD radiation to all 



30 



generator makes use of the Color Dipole Model 

orders. Both use string fragmentation |^ to generate the hadronic final state. 

A complete Monte Carlo simulation of the HI detector response is performed for NC 
and CC DIS processes. For NC DIS, a sample corresponding to Q^ > 1000 GeV^ and 
which amounts to ~ 25 times the integrated luminosity of the data was produced using 
ARIADNE. A similar sample was produced using LEPTO. To reduce further the Monte 
Carlo statistical uncertainty at high Q^, a simulation corresponding to ~ 100 times the 
integrated luminosity of the data was performed using ARIADNE with the additional 



requirement Ee > 50 GeV on the energy of the scattered positron. The same event 
generator was used for the simulation of CC DIS events with Q^ > 100 GeV^ to give a 
sample also ~ 100 times the integrated luminosity of the data. For the comparison with 
experimental data in the following sections, the ARIADNE model for the hadronic final 
state is used unless explicitly otherwise stated. 

The following experimental errors are propagated as systematic errors on the mean 
standard NC DIS expectation: 

• the uncertainty on the integrated luminosity (±2.3%); 

• the uncertainty on the absolute calibration of the calorimeters for electromagnetic 
energies (±3%); 

• the uncertainty on the absolute calibration of the calorimeters for hadronic showers 
(±4%); for NC DIS, this error only enters via the constraints on Premiss and on 

nE-Pz). 

In addition, a 7% "theoretical" uncertainty on the predicted NC DIS cross-section origi- 
nates from contributions of: 

• the uncertainty of 5% on parton density distributions extracted from "QCD fits" ; 
this is partly due to the experimental errors on the input data (in particular in 
the structure function measurements in the high x range from the BCDMS exper- 
iment 0]) and is partly linked to the assumptions for the shapes of the parton 
distributions at the Qq at which the perturbative QCD evolution is started; this 
uncertainty is compatible with what can be inferred from a comparison of the cross- 
section predictions, calculated with HECTOR [^ , using the recent sets of next-to- 



leading order parametrizations of the MRS ||33|, CTEQ ^^, and GRV ||3^ groups 



regulated in the MS scheme |^ 



the value of the strong coupling constant a^, which leads to an uncertainty of 4%; 
this was inferred by comparing the CTEQ (Al) to (A5) |^ parametrizations with 
as{Mz) ranging from 0.110 to 0.122; a similar value can be inferred by a comparison 
of the MRS (Rl) and (R2) ||3^ parametrizations; 

the higher order QED corrections imply a 2% uncertainty in the y range considered 
in this analysis; this was estimated using HECTOR |^ which makes possible two 



approaches: QED radiative corrections can be calculated PSI] in the leading logarith- 
mic approximation at 0{aem) including 0{al^) corrections in the next-to-leading 
logarithmic approximation as well as soft photon exponentiation; alternatively |3^ , 
complete QED and electroweak corrections are calculated at 0{aem.)- 

All analyses described in the following sections have been repeated with an independent 
shift of the central values by ±1 standard deviation of each of the experimental and 
theoretical sources of errors. The overall systematic error of the standard DIS model 
prediction is determined as the quadratic sum of the resulting errors and of the statistical 
error on the Monte Carlo simulation. 



5 Background Sources 

The contributions from background processes which could give rise to events with true or 
misidentified isolated positrons at high Et or to events with a large PT,miss were evaluated 
using the full simulation and reconstruction program chain. The list of processes is given 
in Table |l[ 

The production of real electroweak vector bosons Z^ and W^ was modelled using the 
EPVEC event generator [^]. This includes processes where the partonic structure of a 



photon is resolved. It should be noted that inelastic single Z° and W^ production is not 
contained in the standard DIS model described in section ^ and thus it is treated here 
as a background source. Especially contributions from radiation of a Z° or W^ from a 
quark can give rise to forward positrons and thus can mimic NC DIS events at very high 

Direct and resolved photoproduction processes were modelled using the PYTHIA gen- 



erator ||4T|]. It is based on leading order QCD matrix elements and includes initial and 



final state parton showers calculated in the leading logarithm approximation, and string 



fragmentation ||31|. To enhance the generated yield of events in the hard scattering re- 
gion, the transverse momenta pT in the hard sub-processes are required to exceed 10 GeV 
for light quark flavour production, 3 GeV for prompt photon processes, and 8 GeV for 
heavy quark flavour production. These cuts do not affect the production of jets with 
Pt > 25 GeV (or of jets containing hard fragments with J2e,-y Pt > 20 GeV), or the pro- 
duction of prompt photons with Pt > 10 GeV. The renormalization and factorization 
scales are both set to P^- The GRV (G) leading order parton densities in the photon are 
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Contributions from two-photon (77) processes, where one 7 originates from the proton, 
were also considered. Electron-positron pair production e~^p — >■ e"*" -|- e~^e~ + X, where 
at least one of the leptons has high Pt, was simulated using the LPAIR generator |43| . 
LPAIR includes the relevant Bethe-Heitler diagrams to lowest order in aem (i-e. order ctg^) 
but does not include the photon bremsstrahlung diagrams where the photon converts into 
a lepton pair or any diagrams in which real or virtual photons are replaced by Z° bosons. 



The contribution from qq production in 77 collisions was simulated using PYTHIA [41 



by folding in appropriate 7 fluxes |^^ , and it did not include diagrams involving massive 
electroweak bosons. The diagrams for both final states {e~^e~ and qq) which involve the 
formation of a Z^ vector boson are included in the calculations of electroweak vector 
boson production using the EPVEC generator. 

When the high Q^ DIS selection criteria which are described in section ^ are applied 
to the above background sources, the total remaining background contamination is found 
to be below 1% (95% CL). Nevertheless, in order to ensure that the background contam- 
ination is negligible everywhere in the kinematic plane compared with the uncertainty 
on the standard DIS expectation, additional cuts are imposed. Each cut is specifically 
designed to eliminate a given background source. 

Background contamination in the NC DIS selection is suppressed by the following: 

1. processes leading to multi-lepton final states are eliminated by requiring at least one 
reconstructed jet with -Erjei > 7 GeV in the polar angular range 10° < 9jet < 145°; 



Partonic Process 
(example) 



Generator 
Model 

[refs.] 



Simulated 

Luminosity 

( Pb-^) 



Upper Limits (95% CL) 

on Events per 

14.19 pb-i 



Single W and Z Boson Production 



e + q 



W 



e + X 
e + V 
jet + jet 



m\ 



8600 
1400 



0.005 
0.16 



e + q 



Z 



e + X 

r+ + r~ 
jet + jet 



m 



42000 
59000 
3600 



0.003 
0.001 

0.18 



High Pt jets photoproduction (direct and resolved) 



7 + g 



q + q' 



1 + 9, 1 + 9 

(hard jet) 
(hard fragmentation) 



41 



500 
500 



0.38 
0.31 



Heavy Flavour Production (direct and resolved) 



1 + 9 
1 + 9 



c + c 
b + b 



41 



600 
600 



0.07 
0.07 



Prompt Photon Production (direct and resolved) 



1+q^q+l 



41 



500 



0.09 



Two-photon processes 



7 + 7 — i> e+ + e 
7 + 7 — >q + q 



43 



4500 
17000 



0.04 
0.002 



Table 1: Upper limits (95% CL) on the expected number of events from background processes 
which survive the event selection for NC DIS candidates and the full set of background rejection 
cuts; for the photoproduction of high Pt jets involving light quark flavours, two alternative 
selections are made which require either a hard jet or jets containing hard leptons, photons or 



vrO's. 



the jet is specified using a cone algorithm with a radius JAriJ^f + A^^g^ = 1; at 
least 1% of the jet energy should be deposited in the hadronic section of the LAr 
calorimeter; 

2. QED Compton events are suppressed by rejecting events in which a second positron 
(or photon) cluster with i?Te(2)/-Sr,e > 0.9 is found back-to-back in azimuth within 
A0e,2 > 160°; 

3. background from single Z° boson production followed hj a Z^ —>■ e~^e~ decay is 
suppressed by rejecting events in which the invariant mass of two electron/positron 
clusters lies within 5 GeV of the Z° boson mass; 

4. background from single Z^ boson production followed hj a Z^ ^ t^t^ decay is 
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eliminated by the requirement J2{E — Pz) > 0.2 x Px^miss + 43 GeV, which rejects 
events having both abnormally high PT,miss and abnormally low !](-£' — Pz)'-, 

5. background from isolated prompt photons in photoproduction processes is elimi- 
nated by requiring that there be at least one charged track within the positron 
isolation cone; 

6. background from photoproduction and low Q^ NC DIS are removed by the require- 
ment that there be less than 5 GeV in the backward calorimeters; moreover there 
should be no energy E^-tag in the e-tagger, unless | E^-tag + E^_tag ~ -E'e I < 5 GeV, 
where E^^tag is the energy in the 7-tagger, in order to save NC DIS events which 
coincide randomly with Bethe-Heitler bremsstrahlung; 

7. background from e+e~ pair production in 7-7 collisions is strongly suppressed by 
requiring that there be in addition to the "scattered" positron at most one elec- 
tron/positron candidate with Et > 4 GeV. 

The cuts (2), (3), (4), (5) and (7) each imply individually an efficiency loss for NC 
DIS processes below 1% while the corresponding losses for cuts (1) and (6) are below 
2%. Globally, the complete set of background rejection cuts leads to a 5% efficiency loss 
for NC DIS. It should be noted that the background originating from single Z^ or W^ 
production followed by hadronic decay of the boson, and leading to multi-jet final states, 
could be considerably reduced by requiring that at most one jet carries Exjet > 40%i?T,e- 
This requirement has not been imposed in this analysis in order to stay as close as possible 
to an inclusive DIS measurement. Applying all selection cuts and background rejection 
cuts, the remaining background contamination is below 0.1% (95% CL) for the NC DIS 
sample. Hence, the background is negligible when compared to the estimated errors on 
the standard NC DIS expectation. This statement is also valid in a restricted kinematic 
range of Q^ > 10000 GeV^ where less than 0.1 (95% CL) background events are expected. 
In most cases, the estimation of upper limits for the various background sources in Table |l] 
are limited by available Monte Carlo statistics. 

To deal with specific background sources to CC DIS, it is required that there should be 
no isolated positron with Et^c > 10 GeV within 10° < 9^ < 145°. This rejects in particular 
background from single W boson production followed hj a W ^ eu decay. This causes 
negligible efficiency losses for the CC DIS selection, while the remaining contamination 
from single W production is below 0.02% (95% CL) within the analysis cuts. 

The purity of the event selection is cross-checked by two methods: 

1. for the NC selection, the positron isolation criteria were relaxed and the positron 
cluster to charged track link was not applied; no evidence for background was found; 

2. all events from the NC and CC selection were scanned visually; in the selected NC 
sample, no contamination or overlap with cosmic ray showers or beam halo events 
was seen; there was no evidence for wrongly identified positron clusters; the selected 
CC sample showed no sign of contamination due to cosmic ray showers; no muon 
candidate was found in any of the selected events. 
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6 Analysis 

6.1 General Properties of NC and CC Events 

The measured data are compared with the expectation from the standard DIS modeL Here 
and in all that follows, the expected yields of events are calculated using Monte Carlo 
simulation of the HI detector response, including all detector and selection efficiencies, 
and the full HI event reconstruction program. All comparisons with measurement are 
made with a normalization specified by the total cross-section within acceptance and the 
measured integrated luminosity. 

We find in the experimental data 443 (460) NC DIS candidates satisfying the selection 
criteria and Ql > 2500 GeV^ (Qla > 2500 GeV^), which is in good agreement with the 
expectation of 427 ±38 (442 ±40) from NC DIS processes. Within the kinematic range in 
6e, ET,e, He and Ql considered, the ratio of this expected number of events to the number 
generated by the Monte Carlo simulation is about 0.8 due to experimental acceptance 
and event selection. The sources of errors for the NC DIS expectation are discussed 
above. For the CC DIS candidates, 31 events are found satisfying the above requirements. 



c 
!5 

« 10' 

c 
<u 
> 
<u 

10 



10 




100 200 300 
E rceV) 



c 
!5 

g lO'' 

0) 



10 



• HI data 
-NCDIS 



(b) 




25 50 75 100 
0, (deg) 



• HI data (c) 
-NC DIS 




C 35000 






'• 


O 30000 
Q^ 25000 


; (d) 


• • 


• 


20000 


- 


• 
• 


- 


15000 


r 


*•• 


-_ 


10000 


J 


? 


-_ 


5000 


/ 




-_ 


n 









10 20 30 40 
PT.„,iss (GeV) 



10000 20000 30000 
Q^ (GeV^) 



Figure 1: For NC DIS candidates, the distributions of (a) the measured positron energy E^, 
(b) the polar angle 9e, (c) the missing transverse momentum PT,miss, and (d) the correlation of 
QIq a^d Qg, (symbols); in (a), (b) and (c) the expectations from standard NC DIS assuming 
the integrated luminosity of the measurement are shown as superimposed histograms. 
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which agrees well with the expectation of 34.2 ± 5.8 from CC DIS processes. Here the 
corresponding ratio of accepted to generated events is again about 0.8 in the PT,miss and 
yh range considered. 

Fig. |1] shows, for NC DIS candidates, the distributions of energy and polar angle of 
the scattered positron, the missing transverse momentum, and the correlation of the 4- 
momentum transfer squared Q^ calculated with the e-method (Ql) and the 2a-method 
(QD- Superimposed as histograms are the expectations from the NC DIS simulation. 
The distributions are well described by the simulation except for the largest values of the 
measured positron energy Ef.. This will be quantified below in terms of Q^ to which E^ 
is closely related at large values. A good correlation is seen between the measurement of 
Q^ by the electron and the 2a-method. 

Fig. ^ shows the measured and expected differences between various mass reconstruc- 
tion methods. The e-method, which provides good resolution for all kinematic quantities 
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Figure 2: Distributions for the NC DIS sample of the differences between the mass (M = y/xs) 
calculated using the e-method (Me) and (a) using the 2a-method (M2q), (b) using the /i-method 
(Mh), and (c) using the invariant mass of the positron-jet system [Me+jet)'-, superimposed on the 
data points (• symbols) are histograms of the expectation from the Monte Carlo simulation; (d) 
difference between Me and the generated mass M^en when the latter is calculated either from 
the final state positron using simulated NC DIS events (histogram), or taken from the invariant 
mass of the positron-quark sub-system in a specific eq resonance model (dotted histogram). 
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Figure 3: Selected NC DIS candidate events in the M^ - ye plane; three contours of fixed Q^ 
are shown. 

at tje > 0.1, is taken here as a reference. Good agreement is found between the mass 
measurement Me based on the e-method and M2a based on the 2a-method (Fig. ^ja). 
The difference Mg — M2a is found to be described acceptably by the NC DIS simulation. 
The measurement Mg is compared to Mh calculated from the hadronic energy flow in 
Fig. 1^. Also here good agreement is found between the two methods and the shape of 
the Me — Mfi distribution is well described by the simulation. This confirms the validity 
of the /i-method which is the only one available to reconstruct the kinematics of CC DIS 
events. A comparison of Me with the invariant mass Me+jet calculated from the energies 
and angles of the final state positron and highest Et jet is shown in Fig. ^. Good agree- 
ment is found between the two mass methods which reflects the fact that the measurement 
is dominated by hard e + par ton — > e+parton scattering leading most often to a hadronic 
final state with a single high transverse energy jet. Also the difference Me — Me+jet is well 
described by the simulation. 

Fig. ^ shows the expected experimental contribution to the resolution on Me ob- 
tained by comparing for the NC DIS simulation the reconstructed and generated positron 
momentum 4-vectors. The mean value of the distribution of the reconstructed Me is 
within ±1% of the true value of the Lorentz invariant M = ^/xs calculated from the gen- 
erated final state positron for the range of y considered. With respect to this generated 
positron, values of Me are measured with an RMS resolution of aRMS — 7 GeV and a 
Gaussian peak resolution of (Tcauss — 3.0 GeV. Also shown as a dotted histogram is the 
expected resolution of the e-method for the invariant mass of the positron-parton system 
in a specific model (see section |^). It should be recalled that the Me resolution severely 
degrades at low i/e and the events in Fig. ^ are restricted to i/e > 0.1. A shift equal to the 
experimental systematic error on the calorimeter energy scale for positrons of 3% would 
lead to a y dependent shift for the mass Mg of SMg/Me = 0.03/{2ye). 

Fig. ^ shows the two dimensional distribution of j/e against Me. The cut-off in event 
density for low ye or Me arises because of the requirement Ql > 2500 GeV^. Also shown 
are the Ql = 10000 GeV^ and Ql = 15000 GeV^ contours. 

Fig. ^ and |b show the projected Me and ye distributions of the selected events at 
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> 



"low" g2 (2500 GeV^ < Ql < 15000 GeV^) and Fig. |c and |d at "high" Q^ (q 
15000 GeV^). The distributions of the data are well reproduced by standard DIS predic- 
tions in the low Q^ range. At high Q^ the data exceed the NC DIS expectation, especially 
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Figure 4: Distributions of Mg and ye for selected NC DIS candidate events, (a) and (b) for 
2500 GeV^ < Ql < 15000 GeV^ (c) and (d) for Ql > 15000 GeV^; distribution of ye (e) for 
100 < Me < 180 GeV and (f) for Me > 180 GeV; superimposed on the data points (• symbols) 
are histograms of the expectation of standard NC DIS. 

at Me ~ 200 GeV. Fig. ^ and ^f show the i/e distributions for Mg values below and above 
180 GeV. At high Mg, a difference between experiment and expectation is apparent at 
large ye- 

6.2 Q^ Dependence 

6.2.1 Neutral Current Sample 

Fig. 1^ shows for the NC selection the measured Ql distribution in comparison with the 
expectation from standard NC DIS processes. Also shown is the ratio of the Ql distri- 
bution to the NC DIS expectation. Very similar results are obtained for this ratio for 
minimum values oiymin in the range 0.1 to 0.5. The errors resulting from the convolution 
of the statistical error of the Monte Carlo sample and the systematic errors are correlated 
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Figure 5: (a) Qg distribution of the selected NC DIS candidate events for the data (• 

symbols) and for standard NC DIS expectation (histogram); the arrows indicate the Qg values 
for data entries with Ql > 15000 GeV^ and Mg > 180 GeV; (b) ratio of the observed and 
expected number of events as a function of Q^; the lines above and below unity specify the 
iblcr levels determined using the combination of statistical and systematic errors of the NC DIS 
expectation. 

for different Ql bins and are indicated in Fig. ^. They are shown in this figure (and 
in subsequent figures ^ and |^) as lines above and below unity joining the ±lo" errors 
evaluated at the center of each bin. These errors are dominated by the uncertainty in 
the electromagnetic energy scale of the calorimeter and vary between 8.5% at low Q^ and 
30% at the highest values of Ql- The NC DIS expectation agrees well with the data for 
Ql < 15000 GeV^ while at larger Ql the number of events is in excess of the NC DIS 
expectation. 

To quantify this difference, the numbers of observed and expected events with Ql 
above various Qlnin values are given in Table ||. Also given in Table ^ are the Poisson 
probabilities V{N > Nobs) that in a random set of experiments the number of NC DIS 
events N fluctuates to values equal to or larger than the observed number of events Nobs- 
The systematic error 6b on the mean number of expected events b is taken into account 



16 



QLn( GeV^) 


2500 


5000 


10000 


15000 


20000 


30000 


Nobs 


443 


122 


20 


12 


5 


2 


Ndis 


426.7 
±38.4 


116.2 
±11.6 


18.3 

±2.4 


4.71 
±0.76 


1.32 
±0.27 


0.23 
±0.05 


V{N > Nobs) 


0.35 


0.35 


0.39 


6 X 10"^ 


1.4 X 10-2 


2.3 X 10-2 



Table 2: Number of observed {Nof,s) and expected NC DIS (Ndis) events with 4-momentum 
transfer squared Q^ above given thresholds {Qlnin)'-> '^[N > Nobs) is the probabihty that the 
number of NC DIS events fluctuates to values equal to or larger than Nof,g in a random set of 
experiments. 



by using the convolution : 

V{N > Nobs) = / dxGix; b, 6b) J2 ^(^5 ^) 



k=Nobs 



where the following notations have been introduced : 



• p{k; x) is the Poisson probability to observe k events when the number of expected 
events is x, i.e. p{k;x) = e~^x^ /k\; 

• G{x; b, 5b) is the probability density function for the NC DIS expectation x, namely 
a Gaussian of mean value b and width 5b. 

The resulting probabilities are about 1% at the largest Q.^j„ values. For Q^ > 15000 GeV^, 
the number of observed events is Nobs = 12 for an expectation of 4.71 ±0.76 corresponding 
to a probability V{N > Nobs) of 6 x 10-^ 

6.2.2 Charged Current Sample 

Fig. 1^ shows for the CC selection the measured Q\ distribution in comparison with the 
standard DIS CC expectation. The systematic errors are relatively large and dominated 



QL„( GeV2) 
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31 


24 


10 
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3 


Ndis 


34.2 

±5.8 


21.1 

±4.2 


5.07 
±1.88 


1.77 
±0.87 


0.74 
±0.39 


V{N > Nobs) 


0.64 


0.31 


7 X 10^2 


0.14 


5.4 X 10-2 



Table 3: Number of observed (Nobs) and expected CC DIS events (Ndis) with 4-momentum 
transfer squared Q^ above given Q^j„ values; V{N > Nobs) is the probability that the number of 
CC DIS events fluctuates to values equal to or larger than Nobs in a random set of experiments. 
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Figure 6: (a) Ql distribution of the selected CC DIS candidate events for the data (• symbols) 
and for standard CC DIS expectation (histogram); (b) ratio of the observed and expected number 
of events as a function of Q|; the lines above and below unity specify the it la levels determined 
using the combination of statistical and systematic errors of the CC DIS expectation. 

by the uncertainty on the hadronic energy scale of the calorimeter. They vary between 12% 
at Ql ^ 5000 GeV^ and 60% at Ql ~ 30000 GeVl Within errors, the distribution of the 
measured events is reproduced both in shape and in absolute normalization. Table ^ gives 
the number of observed and expected CC events, as well as the Poisson probability V{N > 
Nobs) as defined in section |6.2.1| , for increasing values of Q^m- I^ the kinematic region 
Q\ > 15000 GeV^, there are Nobs = 4 observed events compared with an expectation 
of 1.77 ± 0.87 from standard CC DIS. This small difference corresponds to a Poisson 
probability V{N > Nobs) of 14%, including systematic errors, that the CC DIS signal 
fluctuates to values equal to or larger than Nobs in a random set of experiments. 

6.3 Mass Dependence as a Function of y 

Fig. ^ shows the measured and expected Me distribution and Fig. |^ the ratio of the 
measured Me distribution to NC DIS expectation for a minimum i/e value of ymin = 0.2. 
Similar distributions are shown in Fig. |^c and 0d for ymin = 0.4. An excess of events over 
the NC DIS expectation at the highest mass (~ 200 GeV) is seen, which becomes more 
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Figure 7: (a) Distribution of Me for the observed (• symbols) NC DIS candidates with 
ye > 0.2; the expectation from standard NC DIS is shown as the superimposed histogram; (b) 
ratio of the observed and expected numbers of events as a function of Me for ye > 0.2; the hues 
above and below unity specify the itlo" levels of uncertainty for the standard NC DIS expectation 
using the combination of statistical and systematic errors; (c) as in (a) but now with ye > 0.4; 
(d) as in (b) but now with ye > 0.4. 

visible with the larger ymin cut. 

To quantify the difference in the Me distribution between the data and the expectation, 
the number of events with Me and y^ above various sets of minimal values Mmin and ymin 
are given in Table ^. Also given are the Poisson probabilities V{N > Nobs) that the 
standard NC DIS signal N fluctuates to values equal to or larger than the number of 
observed events in a random set of experiments. The errors on the NC expectation as 
well as the probabilities take into account all systematic errors described in section ^. It 
is noted that the values obtained for the various sets of Mmm and ymin cuts are correlated. 
Good agreement is observed for low Mmin and ymin, whereas at high Mmin and ymin the 
probabilities V{N > Nobs) become smaller. 

The range of Me values for which the most significant excess over NC DIS expectation 
exists is investigated in detail by considering "windows" of various total widths AMe. 
The central values Me of the mass windows are varied in steps of 1 GeV between 80 and 
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Ndis 

V 


8 

4.65 ±1.68 

0.15 


6 
2.19 ±0.48 
3.1 X 10-2 


6 
1.07 ±0.23 
1.2 X 10-3 


4 
0.68 ±0.13 
5.9 X 10-3 


195 Nobs 

Ndis 

V 


6 

3.58 ±1.43 

0.19 


5 
1.73 ±0.40 
3.8 X 10-2 


5 
0.85 ±0.19 
2.4 X 10-3 


4 
0.53 ±0.11 
2.5 X 10-3 


200 Nobs 

Ndis 
V 


4 

2.72 ±1.10 

0.30 


3 

1.32 ±0.33 

0.15 


3 

0.63 ±0.15 
2.8 X 10-2 


3 

0.39 ±0.09 
8.2 X 10-3 



Table 4: Numbers of observed {Nobs) NC DIS candidates and expected (Ndis) standard NC 
DIS events which satisfy the NC DIS selection for different minimal requirements ye > ymin and 
Afg ^ Mmin', also given is the probability V = V{N > Nobs) that the number of DIS events 
fluctuates to values equal to or larger than Nobs in a random set of experiments. 



250 GeV and the numbers of observed and expected events are determined for different 
Umin values. For each mass window, the Poisson probability V{N > Nobs) is determined, 
including the propagation of all systematic errors described above. The probability for 
each value of Me reflects the level of agreement between data and expectation for an a 
priori choice of Mg. To avoid the inevitable discontinuities when values of V are displayed 
due to the small number of observed events, a probability V is calculated averaging over 
3 steps of 1 GeV around the given central Me value. The resulting Poisson probabilities 
as a function of the central value Mg are shown in Fig. || for a representative choice of 
AMe and ymin settings. This 3 GeV range corresponds to the Me spread due to the flnite 
resolution on the measurement of the scattered positron energy (see Fig. ^jd). 

For low ymim the probabilities fluctuate statistically as expected. At masses above ~ 
220 GeV no data event is observed within the mass windows, and the probabilities V{N > 
Nobs) approach unity by construction. There remain ~ 0.3 NC DIS events expected in 



20 



o ' 



Al 
iQ. 



10 



10 



10 




100 150 200 250 
M, (GeV) 



■° 1 

o ' 



Al 

^10 



10 



10 



10 



-1 


\r\. 


/^ 


-2 


>> g 


1 


-3 
-4 


lAM^ = 25 GeV 

:ye>0-3 


HI 









100 150 200 250 
M, (GeV) 



■g 1 



Al 

|10 



-1 



10 



-2 



10 



-3 



10 



-4 



r(c) \ 

lAM^ = 25 GeV 

:y^>0.4 


i ! 

[ ] 

/h1: 



100 150 200 250 
M^ (GeV) 



■S 1 



Al 

|10 

iQ. 



-1 



10 



-2 



10 



-3 



10 



-4 



Kaa/V^ 


/ : 


"\J 1 


J 


r(d) y 


- 


^AM^ = 40 GeV U 

: Ye > 0.4 V 


HI 







100 150 200 250 
M^ (GeV) 



Figure 8: Probability that the standard NC DIS signal fluctuates to values equal to or larger 
than the number of observed events in a random set of experiments for various ymin values. The 
probability is calculated for different sliding mass windows of widths AMg and averaged over 
3 steps of 1 GeV. At the highest masses, where no data events are observed, the probabilities 
approach unity by construction. 

this high mass region for ymin = 0.4. For ymin values of 0.3 and 0.4 and for a window 
width varying between 25 and 40 GeV, the lowest probabilities of 10~^ to 2.6 x lO"'^ are 
observed for a central mass value M^ ~ 200 GeV. For ymin ~ 0.4, the probabilities are 
below 8 X lO""^ for AM^. = 25 GeV and central Mg values chosen a priori in the range 
~ 198 GeV to ~ 208 GeV. This statement holds with or without averaging V over up to 
5 steps of 1 GeV. It should be noted that the systematic error on the energy calibration, 
which implies an uncertainty on the number of expected NC DIS events and is taken 
into account in the statistical calculations, gives rise to an uncertainty of approximatively 
±5 GeV for the position of the minima in Fig. ^j. The results obtained for a fixed central 
mass value of 200 GeV with ymin = 0.4 and for various choices of AMg, are given in 
Table |^. 

In order to estimate the likeliness that in a random experiment a value of V smaller 
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AMe ( GeV) 


20 


25 


30 


40 


Nobs 


5 


7 


7 


7 


Ndis 


0.63 
±0.13 


0.95 
±0.18 


1.10 
±0.19 


1.57 
±0.28 


V{N > Nobs) 


5.0 X 10"^ 


2.6 X 10"^ 


2.5 X 10"^ 


1.6 X 10"^ 



Table 5: Number of observed NC DIS candidates (Nobs) and expected standard NC DIS 
events (N^js) for different mass windows of total width AMg, an a priori choice for the central 
mass value of 200 GeV and for ye > 0.4; also given are the probabilities V{N > Nohs) as defined 
for Fig. 0. 



than the observed probabilities is obtained anywhere in the mass range from 80 to 
250 GeV, a large number of Monte Carlo experiments were performed. For each of these 
experiments, events were randomly chosen according to the NC DIS expectation. The 
mean number of events in these experiments was conservatively taken to be equal to the 
number of observed events in the corresponding M^ and i/e range. Applying the same slid- 
ing mass procedure as used above for the real experiment and comparing with NC DIS 
expectation including error propagation, less than 1% of all Monte Carlo experiments 
yielded a minimum value V below those obtained in Fig. ^jc and ^. 

6.4 The Events at Very High Q^, M and y 

We find 7 events in the kinematic region M > 180 GeV and y > OA compared with an 
expectation from NC DIS of 1.83 ± 0.33 (1.75 ± 0.32) using the e-method (2a-method) 
These events are largely responsible for the excess of observed events over standard DIS 
expectation in the Q^, M and y distributions in Fig. ^ ^ and 0. They are listed in Table ^ 
and their properties are discussed here in more detail. One of these events is displayed in 
Fig. |. 

For most of these events, the positron is well contained within the fiducial volume of 
the LAr calorimeter. The corrections for material in front of or in between calorimeter 
modules, which are applied during the reconstruction of the energy from showers in the 
calorimeter, were found to be less than 1% (9%) except for the electromagnetic (hadronic) 
showers of events 3 and 6 (6). The electromagnetic showers initiated by the positron of 
events 3 and 6 are corrected by 23% and 5% respectively because they develop near a 
projective azimuthal edge of a calorimeter module. The corrections to the energy of the 
hadronic showers of event 6 is 29% also largely due to shower leakage into dead material 
in between calorimeter modules. The agreement between the different estimators of the 
kinematic variables in Table ^ demonstrates that the corrections are well understood on 
average. In particular, the two estimators Me and M2a are seen to agree well with each 
other and also with the invariant mass Me+jet calculated from the positron and the highest 
Et jet in the final state. 

In each event, the identified positron cluster is geometrically linked to at least one 
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Figure 9: An event of the NC DIS sample at Me > 180 GeV and ye > 0.4. 

high momentum charged track. In one case (event 4), multiple tracks or track segments 
are associated spatially with the positron cluster as expected from Monte Carlo studies 
in 15 ± 5% of events with similar kinematics, given the inactive material between the 
calorimeter and the ep interaction vertex. The 7 events are not confined to any particular 
azimuthal region of the detector. One event (5) is found to have 2.9 GeV deposited in 
the photon detector used to determine luminosity from Bethe-Heitler interactions. The 
probability that one event out of seven of the high Q^, high y sample be in random 
coincidence with an elastic Bethe-Heitler interaction with one photon detected is estimated 
to be ~ 35%. None of the 7 events have high momentum muons which penetrate the 
calorimeter and instrumented iron filter. For each of the 7 events, only one jet is found 
with ETjet > 15 GeV with the cone algorithm and this jet carries more than 90% of the 
transverse hadronic energy flow. 

All the above indicate that these 7 events are well measured NC DIS-like candidates. 
From Table |^, the weighted average mass < Mg > measured for these candidates is 
200.8 ±2.2 GeV. This value does not include the systematic uncertainty of 3% associated 
with the absolute energy scale of the electromagnetic calorimeter. This uncertainty leads 
to a correlated error for all events which depends on i/e like SMe/Me = 0.03/ {2ye), Sye = 
0.03(|/e - 1) and SQl/Ql = 0.03. 

For the CC DIS sample, in the kinematic region Q\ > 15000 GeV^, which includes the 
region in which an excess is observed in NC DIS, there are Nobs = 4 events in agreement 
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Run # 


Me 


M2a 


Me+jet 


Ve 


y2a 


Ql 


QL 




Event # 


(GeV) 


(GeV) 


(GeV) 






( GeV^) 


( GeV^) 


1 


85528 
71329 


196 

±5 


198 


195 


.439 

±.014 


.434 


16950 
±360 


17100 


2 


87050 
8062 


208 

±4 


200 


213 


.563 

±.012 


.582 


24350 
±430 


23320 


3 


88999 
106218 


188 
±12 


185 


184 


.566 
±.032 


.573 


19950 
±1400 


19640 


4 


119314 
18272 


198 

±2 


199 


197 


.790 
±.008 


.787 


30870 
±530 


31320 


5 


122145 
69506 


211 

±4 


227 


211 


.562 

±.012 


.526 


25030 

±440 


27100 


6 


155768 
67216 


192 
±6 


190 


215 


.440 
±.016 


.443 


16130 
±400 


16050 


7 


158577 
100110 


200 

±2 


213 


210 


.783 
±.008 


.762 


31420 
±540 


34450 



Table 6: NC DIS candidates satisfying the kinematic requirement M > 180 GeV and y > 0.4 
when M and y are calculated with either the e-method or the 2a-method; the errors given for 
the e-method take into account the energy and angular resolution for the positron, and a 30% 
uncertainty assigned to dead material corrections; the systematic uncertainty of 3% associated 
with the absolute energy scale of the electromagnetic calorimeter (see text) is not included. 



with the expectation of 1.77 ± 0.87 from standard CC DIS. These 4 events are listed in 
Table 0. 

7 Non-Standard Contributions at High Q^ 

Having presented the essential experimental evidence, we now discuss possible contribu- 
tions to an excess of events at high Q^ (high M and y) within and beyond the Standard 
Model. 



7.1 New Contributions in the Standard DIS Model 

The standard DIS expectations given in section ^ take into account a systematic error due 
to a specific choice of the parton density parametrizations. No parton density parametriza- 
tion compatible with existing experimental DIS data was found to lead to an enhancement 
of the NC DIS expectation at high masses (M > 180 GeV) for Q^ > 15000 GeV^ beyond 
the systematic uncertainty quoted in section ^. This is in particular the case for the 
parametrization CTEQ (4HJ) [Q which was designed to cope with the apparent rise of 
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Run # / Event # 


M, 


Vh 


Ql 


1 


85987/ 99058 


213 


.52 


23380 


2 


153720/199055 


187 


.58 


20170 


3 


163852/20613 


157 


.75 


18650 


4 


169851/206239 


192 


.56 


20800 



Table 7: CC DIS candidates satisfying the kinematic requirement Q\ > 15000 GeV ; in the 
M, y and Q^ range spanned by these events, the expected resolutions are 5Mh/Mh ~ 10%, 

associated with 



^yh/Vh ~ 10%, SQf^/Qj^ ~ 28%, not including the systematic uncertainty of ■ 
the absolute hadronic energy scale. 



the di-jet cross-section observed at the Tevatron [^ and leads to an enhancement, com- 
pared to MRS (H), of < 5% in the above cited M and Q^ range. In the same context, the 
parametrization MRS (R2) [Q which assumes a "large" as value {ag^M^) — 0.120) leads 
to an enhancement of the scaling violation and consequently to a slight reduction of the 
differential cross-section at high M relative to MRS (H), also within systematic errors. 



The MRS (J') parametrization [^ attempts to reproduce the magnitude of the rise 
of the di-jet cross-section observed in particular in CDF |4^ but does not fit existing DIS 
data (e.g. BCDMS data at high x and low Q^). It would lead to a cross-section about 
20% larger than the values used for this analysis which would still not explain the excess 
of events observed in the above cited M and Q^ range. 

The statistical significance presented in previous sections for an excess at high Q^ 



and/or high M and y is unchanged when using the LEPTO [Q model for the QCD 
corrections. 

The expectations from the standard DIS model have been given under the conservative 
assumption that the longitudinal structure function Fl is equal to zero. A finite value of 
Fl either originating from QCD or due to Fermi motion in the proton, would decrease 
the cross-section at large values of Q^ and y. 

From the above, and given the agreement of the measurement with NC DIS expectation 
at low M, Q^ and y (see Fig. ^, there appears to be no mechanism within the standard 
DIS model framework which would lead to an enhancement of the cross-section at high 
Q^ or high y as observed. Within the standard DIS model the only explanation of our 
result is therefore a statistical fluctuation. 



7.2 Physics Beyond the Standard Model 

Beyond the Standard Model, the excess at high Q^ and/or high M and y could be ex- 
plained by different mechanisms. 

New particles could be produced as resonances in the positron-parton system. Promi- 
nent examples for new particles with couplings to positron-parton pairs are leptoquarks |^ 



leptogluons |4^ and squarks in R-parity violating versions of Supersymmetry |^. It 
should be emphasized that the reconstruction of the invariant mass of such a resonance 
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is hampered by QED and QCD corrections. An example of such an effect in the form of 
a narrow s-channel resonance is shown in Fig. ^d as a dotted histogram. The LEGO [Q 
generator, which was used to produce such a resonance at fixed mass M^en, incorporates 
initial state QED bremsstrahlung in the collinear approximation and QCD initial and 



final state parton showers and fragmentation [^. It moreover properly takes into ac- 
count the effects of the parton shower masses on the decay kinematics. It can be seen 
that a resonance would be reconstructed with a systematic shift (typically < 5%) towards 
smaller mass values and with a resolution considerably worse than expected from detector 
effects alone. At the present level of significance, neither the RMS width of the measured 
events of 8.5 GeV nor the ye distribution (Fig. |f), which is expected to extend to larger 
He for a resonance than for NC DIS processes, allow conclusions to be drawn concerning 
the possibility of resonance formation. 

A slight excess is seen in the Q^ distribution (Fig. |^, Table 0). The virtual exchange of 
new particles between positrons and partons or possible substructure of fermions inducing 
contact interactions could all lead to an enhancement of the cross-section at high Q 
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8 Summary 

Deep-inelastic scattering (DIS) events have been observed in e^p collisions at very large 
4- momentum transfer squared Q^, and have been compared with the expectations from 
the standard Neutral Current (NC) and Charged Current (CC) DIS model. 

For Q"^ < 15000 GeV^ the distributions of Q"^ oi M = y/xs and y = Q'^/M'^ are well 
reproduced by the expectation of standard DIS. At larger momentum transfer (Q^ > 
15000 GeV^), 12 NC DIS candidate events are observed where 4.71 ± 0.76 are expected 
and 4 CC DIS candidate events are observed where 1.77 ±0.87 are expected. The Poisson 
probability V that the signal from standard DIS fluctuates to a number of events equal 
to or larger than the observed number of events is 6 x 10"^ in the NC case and 0.14 in 
the CC case. 

For the NC candidates, the excess of events is most prominent in a mass window of total 
width 25 GeV centered at an invariant mass M ~ 200 GeV of the positron-parton system. 
This mass is consistently determined using different kinematic reconstruction methods, 
for which only the measured positron, or only the measured angles of the positron and 
the hadronic system, or only the 4-momenta of the positron and the jet with highest Et 
are used. The dominant e + jet + X topological feature of the events is characteristic 
of standard neutral current DIS processes. In a mass window of 25 GeV width with a 
central value of 200 GeV and for y > 0.4, 7 events are observed where 0.95 ± 0.18 events 
are expected. For this and other choices of mass windows and y thresholds, the probability 
to observe an excess as large as the measured one anywhere in the mass range investigated 
is of order 1%. 

No known detector effect can account for an excess at large Q^ or can be associated 
with an excess which occurs preferentially in a restricted range of M. 

Given the existing experimental constraints on parton density distributions at high M 
and lower Q^, and given the agreement of the resulting predictions for e~^p DIS which is 
here reported for Q^ below 15000 GeV^, there is little freedom for an enhancement of the 
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cross-section at higher Q^ by different choices for the partonic structure of the proton, by 
changing the strong couphng constant, or by including higher order corrections. 

A new mechanism would be needed to explain an enhancement of the DIS cross- 
section affecting mostly high Q^ or high y values. Within the standard DIS model the 
only explanation of our result is therefore a statistical fluctuation. 

Whereas an account of the observation by introducing new physics beyond the stan- 
dard model of electroweak and strong forces is kinematically possible, the signature of the 
observed events is identical to DIS and hence a clarification of their nature will have to 
come from the study of kinematic distributions with larger statistics. 
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